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The parity violating longitudinal asymmetry A is calculated for quasielastic elec- 
tron scattering. We use a variety of relativistic mean field models for the response 
of nuclear matter and 12 C at a momentum transfer of q=550 MeV/c. Relativistic 
effects from a reduced nucleon mass, RPA correlations and vacuum polarization can 
all change A by a relatively large amount. These large nuclear structure correc- 
tions could make it impossible to extract strange quark matrix elements or radiative 
corrections to weak axial currents. 
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I. INTRODUCTION 



Parity violating electron scattering probes weak neutral currents and can provide very 
interesting information on strange quark matrix elements in nucleons and nuclei. One is 
interested in the form factors of the the vector (both electric and magnetic) and the axial 
strange quark currents in a nucleon. These three form factors contain new information on 
nucleon structure which can help guide and constrain QCD based models. 

It will take a number of measurements to determine separately all of these form factors. 
Furthermore, there are important complications from radiative corrections which hinder 
the extraction of strange quark matrix elements. Therefore, one anticipates a program of 
several electron scattering experiments involving proton and different nuclear targets. These 
results will then be combined with neutrino scattering in order to extract the most accurate 
strange-quark information. 

Clearly one experiment is elastic electron proton scattering with polarized electrons (e,p) . 
The parity violating analyzing power (the cross section difference between positive and 
negative helicity electrons over the sum) is sensitive to weak neutral currents. Indeed, the 
SAMPLE experiment will place (somewhat crude) limits on the strange quark magnetic 
form factor. 

Neutrino scattering may provide the best information on the axial-vector form factor. 
For example, the BNL experiment measuring neutrino and antineutrino proton scattering 
claimed to determine a nonzero strange quark contribution Note, for experimental 
reasons, most neutrino experiments measure a combination of elastic scattering from free 
protons and quasielastic scattering from nucleons bound in nuclei. In a later paper we will 
examine nuclear structure corrections to quasielastic neutrino scattering and how these could 
affect the extraction of strange quark matrix elements. 

Even if the strange axial current was accurately determined (in neutrino scattering), 
uncertainties in axial currents could still haunt electron scattering experiments. This is 
because of radiative corrections to the axial current which are very different (larger) for 
electrons than for neutrinos [|TJ. For example, a small "wrong" parity admixture could 
exist in the nucleon (perhaps due to weak couplings of the nucleon to its meson cloud). 
The electron could then scatter completely electromagnetically from the small wrong parity 
admixture. These radiative corrections could be large (since they would not involve the 
small weak vector coupling of the electron (see Eq. ( |17c| ) and Table |). Furthermore, some 
radiative corrections are difficult to calculate since they involve unknown strong interaction 
matrix elements. 

One alternative is to measure the radiative corrections by performing another experiment 
(besides elastic electron-proton scattering) which has a different combination of radiative 
corrections and strange quark contributions. One appealing possibility is to measure parity 
violation in quasielastic electron scattering. This is favorable experimentally because it has 
a large cross section and a relatively large asymmetry. Indeed, there has already been one 
pioneering measurement on 9 Be [Q. [However, this experiment was not accurate enough to 
determine strange quark matrix elements]. 

Because quasielastic scattering also involves neutrons one measures a different isospin 
combination than in elastic proton scattering. This could allow one to separate isoscalar 
strange quark contributions from radiative corrections to the isovector axial current. With 
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this in mind, a parity violating quasielastic experiment has been proposed by Donnelly and 
collaborators |5|||. 

However if a nuclear target is involved, one must immediately address two questions: 
(a) do unknown nuclear structure effects complicate the extraction of strange quark matrix 
elements and (b) are the strange quark matrix elements themselves modified in a nucleus? 
We discuss the second (very interesting) question first. 

It is plausible that strange quark currents are different in a nucleus than in a nucleon. 
For example, the old EMC experiment revealed that up and down quark structure functions 
are changed by about 10%. One might expect the strange quark content to increase with 
density. Indeed, in the limit of very high density one expects a transition from nuclear 
to quark matter. Furthermore, it is believed that the ground state of quark matter (in 
weak equilibrium) is strange matter containing a large number of strange quarks. Here the 
large up and down quark Fermi energies overcome the inhibiting effect of the strange quark 
mass. Even at much lower densities, one can speculate that increasing the up and down 
quark Fermi energies would help produce more virtual strange-quark pairs. Thus the strange 
quark content could increase with density. 

In a hadronic picture the same conclusion can be reached. Calculations for neutron 
matter suggest that the fraction of real hyperons increases with density at high density [|7| . 
Furthermore, if a Kaon condensate H were to form at high density, the strange quark content 
would be very large. At lower densities, one could expect the interactions between nucleons 
and the large nucleon Fermi energy will help drive virtual transitions of nucleons to Kaons 
and hyperons. Thus the "Kaon cloud" in a nucleus could be larger than in a nucleon. This 
would increase the strange quark content of nuclear matter. 

At present very little is known about the strange quark content of nucleons much less 
that of a nucleus. Therefore, we will make the minimal assumption that the strange quark 
matrix elements are the same in a nucleus as in a nucleon. We leave the density dependence 
of the strange quark contributions as a very interesting question for further work. 

In this paper we address the first (more pedestrian but nevertheless crucial) question. 
How much do nuclear structure uncertainties complicate the extraction of strange quark ma- 
trix elements from a quasielastic electron scattering measurement? Because the asymmetry 
at backward angles involves a ratio of the parity violating weak response to a very similar 
parity conserving electromagnetic response most nuclear structure effects simply cancel in 
the ratio. We discuss this in detail below. However, the weak vector current has a very 
different isospin character from the electromagnetic current. Therefore, one is sensitive to 
the ratio of isoscalar to isovector responses. A nuclear structure effect that changes the ratio 
of isoscalar to isovector responses will change the asymmetry and destroy ones ability to 
extract strange quark matrix elements. 

Donnely et ai, have argued that this isospin ratio is small and under control (in the 
transverse response) |J. [Note, the transverse response dominates over the longitudinal for 
large electron scattering angles, see below] . Most isoscalar effects are small because of the 
very small isoscalar anomalous moment of the nucleon. This argument is valid even if there 
are big random phase approximations (RPA) or other nuclear structure corrections to the 
bare isoscalar response. The overall isoscalar response is still expected to be small because 
the electromagnetic probe only couples very weakly through the small isoscalar anomalous 
moment. 
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However, there can be sizable contributions of the Dirac F x form factor to the isoscalar 
response. In an earlier paper ||, we argued that these contributions could be significantly 
enhanced by relativistic effects. In the present paper we elaborate on these earlier results 
and examine RPA and vacuum polarization corrections. Furthermore, we extend earlier 
nuclear-matter results and present full finite-nucleus responses. 

Our formalism for parity violating electron scattering is presented in Sec. |H] where the 
asymmetry is expressed as a ratio of several nuclear response functions. Definitions of 
these response functions are presented in Sec. |II A| and results in Sec. |T|. For simplicity 
we perform calculations only for 12 C and nuclear matter at a momentum transfer of q = 
550 MeV/c. Sec. [TV] presents a summary of our findings and our conclusions. We find that 
relativistic and other nuclear structure effects are relatively large and that these may destroy 
the interpretability of a parity violating quasielastic measurement. 



II. FORMALISM 
A. Nuclear Response Functions 

In this section we express the parity violating asymmetry in terms of a number of response 
functions for weak and electromagnetic currents. This formalism has also been presented 
in Refs. PflliOfl . The longitudinal analyzing power, or asymmetry, A is the cross section 
difference for scattering electrons of positive and negative helicity divided by the sum. It 
clearly vanishes if parity is conserved. However, in the standard model, interference terms 
between the weak neutral and the electromagnetic current produce a small but nonzero A. 

In lowest order, the inclusive (polarized) cross section 
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contains the one-photon exchange contribution plus the 7 — Z° interference term (the purely 
weak component is negligible and will not be considered further). The leptonic contribution 
to the cross section is contained in the two tensors 



k^k u + k k v <Ji_iv(k • k ) + %hE^ va pk k 13 



(2a) 
(2b) 



Here h is the electron helicity and gy{g e A ) is the weak- vector (weak-axial) charge of the 
electron (see Table |). The response of the nuclear target is contained in the two hadronic 
tensors and W^. The electromagnetic tensor, defined by 



W^(q,w) = X;<*»l^(q)l*o><*o|^(-q)|*»>%-wO 



(3) 



can be expressed, using current conservation, Lorentz covariance, and parity invariance, in 
terms of two independent response functions 
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where g M = (a;, q) is the four-momentum transfer to the target, and p 1 * and M T are the 
four-momentum and rest mass of the nuclear target respectively. Since the electromagnetic 
tensor is symmetric under the exchange of Lorentz (fx <-> v) indices, the electromagnetic 
contribution to the cross section becomes independent of the electron helicity. Furthermore, 
since weak-interaction effects on the unpolarized cross section are negligible, the parity 
conserving contribution to the cross section is entirely electromagnetic and is given by, 



1 

Cm 



d 2 a h 
dtt'dE' 



[v L S L (u,q) +v T S T (uj,q)) 



(5) 



where a Ar is the Mott cross section, v L and v T are kinematical factors 



Q 1 



2cr 



+ tan 2 (#/2); Q 2 = q 2 - 



(6) 



and S L and S T are the longitudinal and transverse nuclear response functions. These are 
related to W\ and W2 by the following simple relations 



Sr. = Wl 



q 2 q 4 



S T = w 11 + w 2 



-2Wi . 



(7a) 
(7b) 



The parity-violating electron asymmetry, driven entirely by the electromagnetic-weak 
interference term, is contained in the hadronic tensor W^. In order to express the 
electromagnetic-weak tensor in terms of Lorentz-invariant response functions we write it 
in terms of vector- vector and vector-axial contributions, i.e., 



rv -fZ — 'Sz T 7Z ' 



(8) 



where W^(W^" 5 ) arises from the interference of the electromagnetic vector current with 
the vector(axial) component of the weak current, i.e., 



W^(q,w) = X;<*»l^(q)l*o>(*o|^ (-q)l*»>*(w-a; B ) 

n 

^r(q^) = E(*ni^ M (q)|*o)(*o|^ 5 (-q)i^)%- 1 



(9a) 
(9b) 



In particular, the purely vector component has exactly the same Lorentz structure as the 
electromagnetic tensor displayed in Eq.(ffl), i.e., 



W 1 



<r - 



q 2 



+ 



Ml 



(10) 



The vector-axial tensor, on the other hand, must be written in terms of the only available 
pseudotensor that one can construct from p M and g M , namely, 



~ Wa 



(11) 
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The parity-violating part of the cross section can now be obtained by performing the 
contraction of the leptonic tensor l lZ v (Eq.([2bD) with the hadronic tensor Wj*%. This yields, 
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Cm 
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dVl'dE' 



a rh g\ (v l S l (uj, q) + v T S T (u, q)J + g e v v T >S A (uj } q) 



(12) 



where 



v T i = tan(6*/2) 



2cr 



.1/2 



+ tan 2 (#/2) 



and the three additional nuclear response functions are given by 



S T = w? z + w% 

S A ee i (W^ - W^) = W A . 



q_ 

-2W 1} 



(13) 

(14a) 
(14b) 
(14c) 



The factor a r (r = Q 2 /AM 2 ) sets the scale for the magnitude of the parity- violating effects 
and is given by, 



GpM 2 



-3.172 x 10" 



ira 



(15) 



Notice that in the above expression only terms linear in the electron helicity were kept. 
Terms independent of the electron helicity make a negligible contribution to the parity- 
conserving cross section and were neglected. The parity- violating asymmetry, defined as the 
difference of helicity cross sections divided by their sum, is now given by 



_ dcr T — d<J l AAA 

acr T + da. 



where 



A L /a T = g 
A T /a r = g 
A T '/a r = g 



Vr.Sr 



A v L S L + v T S T 

A v L S L + v T S T 

v T 'S A 
v v L S L + v T S T 



(16) 

(17a) 
(17b) 
(17c) 



The hadronic tensors defined above are intimately related to a fundamental many-body 
operator, namely, the current-current correlation function or polarization tensor. This rela- 
tion can be illustrated, for example, in the case of the timelike-timelike component of the 
electromagnetic tensor, 



(18) 
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where the polarization tensor, II 00 , has been introduced in the last line. Consequently, all 
five nuclear response functions can be written in terms of appropriate components of several 
polarization tensors, i.e., 



S T (q,u 
S T (q,u 



= — Imn7q )W ), 



7T 



-±Jro(n£(q,a;) + II£(q,a;)) , 
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= — Im n"(q^) + n-(q,c 



7T 

= — Win£< B >(q,u;) 

7T v 



(19a) 
(19b) 
(19c) 
(19d) 
(19e) 



The advantage of relating all nuclear response functions to a many-body operator, i.e., 
the polarization tensor, is that the responses can be systematically computed using well- 
known many-body techniques (e.g., Feynman diagrams). In particular, this will enable us to 
calculate a nuclear response that goes beyond the simple impulse, or uncorrelated, response 
by including many-body correlation effects. 

The various nuclear response functions can, then, be calculated from imaginary parts 
of polarization insertions given a model for: (a) the hadronic weak and electromagnetic 
currents and (b) the nuclear structure. Hadronic currents are discussed in the next section. 
Subsequently, we calculate the nuclear responses in a number of relativistic mean-field models 
with and without RPA correlations and with various treatments of vacuum polarization. 



B. Electromagnetic and Weak Hadronic Currents 

In this section we describe the electromagnetic and weak hadronic currents that will be 
used in calculating the five nuclear response functions. We start by writing the electromag- 
netic current for the three lightest (u, d, s) quarks 

J£ M = ejurfu + ej^fd + e.&fs . (20) 

For two quark flavors (u and d) one customarily rewrites the electromagnetic current in 
terms of isoscalar and isovector contributions. The analogous expression in the three-flavor 
case is given in terms of two octet (3 and 8) and the singlet SU(3) currents, i.e., 

J^em = e qYq + e 3 q^q + e 8 g 7 ^y q 

= ^yg+^^yg- (21) 

where we have defined (and used) singlet and octet electromagnetic charges 

e = (e„ + e d + e s )/3 , (22a) 
e 3 = (e„ - e d ) , (22b) 
e 8 = (e„ + e d - 2e s )/Vs . (22c) 
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Single-nucleon electromagnetic form factors are now obtained by evaluating the electro- 
magnetic current between nucleon states i.e., 



(N(p's't')\JUT = 0)\N(pst)} = ±(N(p's't')\qr^q\N(pst)) 



U(p's' 



2M 



U(p, s)8 t > 



(N(p's't')\JUT = l)\N(pst)} = (N(p's't')\qY'-^q\N(pst)) 



U{p's' 



2M 



u(j>,s){n 



3)t't 



(23a) 



(23b) 



The weak- vector and weak-axial neutral currents can be written in analogy to Eq. (pip. 
First, however, one rewrites the SU(3)-singlet current in terms of one octet (eight) and a 
pure strange-quark current. This yields, 



° JVC 



X 

2 



(24) 



where the couplings are simple linear combinations of quark weak-vector charges (see Table 
ED, i.e., 



& = (9v ~ 9 d v) , 

e s) = (g u v + 9 d v + 9t 



d ) 



(25a) 
(25b) 
(25c) 



Writing the hadronic current as above, enables one to express nucleon matrix elements of 
the weak- vector current in terms of electromagnetic isoscalar and isovector form factors plus 
a strange-quark contribution (see Eqs. ( [23a| , |23a| )) , 



(N'\J£ C (T = 0)\N) = V3 ev } (N'\JUT = 0)\N) + e^iN'lsYslN) , 
(N'lW = 1)\N) = e^iN'lJUT = 1)\N) , 



(26a) 
(26b) 



where we have assumed that matrix elements of the strange-quark current are equal for 
protons and neutrons. Defining single-nucleon strange (F^ 3 \F^ a) ) and weak- vector form 
factors (Fj 0) , F 2 (0) , F^\ F 2 (1) ) in analogy to Eqs. ( |23a| , [23b|) we arrive at the following useful 
relations between single-nucleon form factors 



j>) = ^F^ + ^'FW 
F 2 (0) = v/3 ^ 0) F 2 (0) + 



1 1 J 1 ) 

2 Sy ± 2 



(27a) 
(27b) 



Finally, one defines, in analogy to Eq. (E4), the weak-axial neutral current in terms of 
octet and strange axial- vector currents and the appropriate weak-axial charges (see Table [11]) 



(28) 



Single-nucleon form factors are obtained (neglecting induced-pseudoscalar contributions) by 
evaluating octet and strange axial-vector currents between nucleon states 
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(N(p's't')\qYl 5 jq\N(pst)) = Gf [u{p> 's')rj 5 U(p, s)] (r,)* 



(N(p's't')\qYl 5 jq\N(pst)) = Gf [U(p' s')Yl 5 U(p, s) 



St't 



(N(p's't')\s^ 5 s\N(pst)} = G { ° } U(p's')YfU(p,s) 5 t , t . 



(29a) 

(29b) 
(29c) 



Isoscalar and isovector single-nucleon axial form factors are now obtained from the following 
relations 



(30) 



In conclusion, the electromagnetic and weak-neutral currents that one is to employ in 
calculating all nuclear response functions have been written in terms of appropriate single- 
nucleon form factors (described in the appendix) and are given by, 



EM 



2M 

f^+zf^^ 



'J ATS 



+ 
+ 



F 1 (1) 7 M + iF^a 



2M 



1 ' 2 2M 



(31a) 

(31b) 
(31c) 



The above equations are a major assumption of our model and prescribe how the hadronic 
currents are taken "off-shell". We also assume that the form factors are unchanged in the 
medium. 



C. Relativistic Mean-Field Responses 

In this section we describe the model used in the calculation of the various weak and 
electromagnetic response functions. Our starting point is a relativistic mean field approxi- 
mation to the Walecka model where nucleons move in strong scalar E s and vector S v mean 



fields [12,13]. In particular, the strong scalar field is responsible for reducing the nucleon 



mass from its free-space value to M* 

M* = M - S s . (32) 

For nuclear matter this is a constant that must be self-consistently determined at every given 
density. In a finite nucleus, however, the effective mass depends on position. It is smaller at 
the center and rises back to its free value M at the surface. 

These mean field models provide perhaps the minimal relativistic description of a nuclear 
target. The strong mean fields are very closely related to the spin-orbit potential of the shell 
model. Furthermore, M* provides a minimal description of the binding energy shift seen in 



the position of the quasielastic peak [ITT 



We illustrate the formalism by presenting electromagnetic, weak vector and weak axial 
response functions for nuclear matter. Most of the formalism, however, remains valid in the 
finite system. The appropriate modifications to be carried out for a finite nucleus and for 
including RPA correlations will be discussed at the end of the section. 
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The basic ingredient in calculating all nuclear response functions is the nucleon Green's 
function (or Feynman propagator). The Green's function describes the propagation of nu- 
cleons in the mean fields. In a self-consistent Hartree approximation the nucleon propagator 
is given by |L3 



G(k) = []k + M* 



1 



k 2 - M* 2 + ie 



+ i^8(k° - E* k )Q(k F - |k|)) = G F + G D 



(33) 



where E£ 



Vk 2 + M* 2 and k F is the Fermi momentum. G F is that part of the propagator 
having the same analytic structure as the free Feynman propagator. The density-dependent 
part of the propagator, G D , corrects G F for the presence of occupied states below the Fermi 
surface. 

The electromagnetic (Eq. (|31aj) ) and weak- vector currents (Eq. ( ftlbD ) have been written 
in terms of Dirac and Pauli (or anomalous) contributions. This implies, that the electromag- 
netic (n"") as well as the weak-vector (n^) polarization tensors will have to be evaluated 
in terms of individual polarizations having Lorentz vector v=7 M and tensor t= a**" vertices, 
i.e., (ignoring isospin labels) 



vt 



*n t 7 



d A k 

d 4 k 
d 4 k 

(2^ 



Ti[G{k)YG{k + q)Y] , 



Tr 
Tr 



G{k)rG{k + q) 



%o vp q P 
2M 



G{k)--^G{k + q) qp 



2M 



2M 



(34a) 
(34b) 
(34c) 



Note that, both, the electromagnetic and the weak- vector responses are driven by the same 
polarizations. The only difference between the two responses arises from the use of either 
electromagnetic or weak-vector single-nucleon form factors. One might expect most nuclear 
structure effects to change Eqs. ( |34a| , |34b| , |34c|) in similar ways. Indeed, within a model, one 
could try to combine the three polarizations and express the responses in terms of the electric 
and magnetic form factors. However, we chose to define our off-shell currents directly from 
Eqs. ( 31a , |31b| ) in terms of F 1 and F 2 form factors. Indeed, we find (in the next section) 
different relativistic effects on IP^ than on n£". Thus, in our model, it is important to break 
the response up into these three pieces (vv, vt, and tt). 

The weak-axial polarization, IP^ 5 , can be similarly written in terms of polarizations 
having either a vector or a tensor vertex, and an axial- vector vertex a=7 M 7 5 , 



iIP ; " 5 

va 



d A k 

d^k 



Tr [G{k)YG{k + q)Yl 



Tr 



G(k) 



2M 



G(k + q )rr 



(35a) 
(35b) 



Note, that the Green's function defined in Eq. [33|, describes the propagation of both 
nucleons and antinucleons. Hence, the above polarization tensors contain, in addition to 
the usual particle-hole excitations, nucleon-antinucleon excitations (vacuum polarization). 
However, vacuum polarization has no imaginary part in the space-like region (q > u) probed 
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in electron scattering. Therefore, vacuum polarization makes no contribution to the uncor- 
rected responses. Vacuum polarization will, however, make an important contribution to 
the correlated response that will be presented in Sec. 

Using the above defined polarizations together with the single-nucleon form factors leads 
to the following expressions for the electromagnetic and electromagnetic-weak polarization 
tensors 

rp; (l) = F^F^U^ + (F«F 2 W + F^F[ l) ) IP;" + F 2 w F 2 (l TL7 w , (36a) 
Wf ] = F®F®n£ w + (F®F® + F^F[ l) ) EF t " w + F 2 w F 2 (l) n t 7 w , (36b) 
n «j*w = pWGWn^sM + F^G^W^ . (36c) 

7Z 1 A va 1 2 A ta V / 

To summarize, the response functions depend on a variety of nucleon electromagnetic and 
weak form factors. These form factors are collected in the Appendix. The nuclear structure, 
on the other hand, is contained in the polarization insertions IE^ , EP t " , ITj^ , IE^" 5 , and 

ta 

The response functions can be calculated in various approximations. The simplest ap- 
proximation is to calculate uncorrelated response functions for nuclear matter. This can be 
done either for M* = M, which gives the free Fermi gas results, or for a smaller value of M*. 
These results were presented in Ref. || and are included below for comparison. Note that in 
these approximations the isoscalar and isovector polarizations are identical. The change in 
the responses arises exclusively from the difference in isoscalar and isovector single-nucleon 
form factors. 

Alternatively, one can calculate the response functions directly in a finite nucleus and 
avoid a local density approximation. The finite-nucleus calculations can be carried out using 
a nonspectral representation of the nucleon Green's function as discussed in Refs. fll4jL5 ,16 



Next, long range correlations can be included by calculating RPA polarization insertions 



Refs. []T4 ,[15 , 16|, 17[| . We employ a residual interaction consisting of only isoscalar sigma 



and omega meson exchange with self-consistent parameters (see Table [IT]] ) chosen to fit 



ground-state properties [|18|| . In particular note that the isovector response is insensitive to 
correlations and thus remains identical to its uncorrelated value. Finally, vacuum polar- 
ization can be included, in addition to particle-hole or core polarization, in calculating the 



RPA response. Further details can be found in Ref. 16 



III. RESULTS 

We now present response functions for the N=Z target 12 C. We expect similar results for 
heavier nuclei. For simplicity we consider only a momentum transfer of q = 550 MeV/ c. Tra- 
ditional nuclear structure uncertainties will be larger for much smaller momentum transfers. 
We expect similar relativistic effects for momentum transfers greater than q = 550 MeV/c, 
but this remains to be investigated. 

First we notice that the parity-violating asymmetry A is dominated by the transverse 
weak vector contribution A T with A T > making perhaps a 20 percent contribution for large 
scattering angles and A L contributing perhaps 10 percent at forward angles. Furthermore, 
A T is itself dominated by the isovector response from the _F 2 (1) form factor contributions (last 
term in Eq. (|36b|) ) . 
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For a large-angle experiment, say at an electron scattering angle of 150 degrees, A L 
is unimportant (as is the contribution from the longitudinal response S L to the parity- 
conserving cross section) because of the kinematic factors v L and v T (v L /v T — > at backward 
angles). Thus, let us focus on A T . We will return to the small A T > contribution at the end. 
Furthermore, we will neglect strange-quark contributions so that, from Eqs. ( |27a| and |27b| ), 
F<® = ^3 and = ^F t w . Thus A T only depends on the ratio of isoscalar S { T 0) 

to isovector S T 1} electromagnetic responses. The only nuclear structure that A T is sensitive 
to is isospin. At large angles then, A T becomes (neglecting the contribution from strange 
quarks) 



(i) 



1 + 



where the isospin ratio is given by, 



Tin 



Tin 



C(0) / cm 

T / T 



(37) 



(38) 



The couplings -\/3 t$ and ££f> (Table ||) have opposite signs and almost equal magnitude. 
Therefore a change in the ratio 1Z T by 0.01 {i.e., TZ T — > TZ T ± 0.01) will change A T by about 
two percent. However, the total transverse response (S^ 0) + S^) is predominantly isovector 
because of the large isovector anomalous moment of the nucleon. Thus TZ T is expected to be 
small (of the order of 5%). However, one must investigate the nuclear structure uncertainties 
in TZ T . Most traditional corrections to the isoscalar response are hidden because of the small 
isoscalar anomalous moment of the nucleon. 

However, there could be large relativistic corrections to the isoscalar transverse response 
A smaller value of M* enhances and thus increases S^K [Note, is dominated by 

rather than by n" because of the small value for F^°\] The smaller value of M* enhances 
a nucleon's velocity and produces a larger convection current (which is half isoscalar). This 
could lead to a large change in A T . 

Below, we present results showing that the uncertainties in the transverse response, and 
hence on A T , are large. For simplicity we do not show results for longitudinal responses. 
However, most of these have already been published in Ref. [TB]. It is already generally 
agreed that nuclear structure uncertainties in the longitudinal response (and hence in A L ) 
are large. 

We present a catalog of transverse responses for ten different calculations. First, we con- 
sider four different finite nucleus calculations. Uncorrelated (or impulse) and RPA responses 
are shown first neglecting vacuum polarization in a so called mean field theory (MFT) and 
then with vacuum polarization in a relativistic Hartree approximation (RHA). For compari- 
son we also show the corresponding four responses calculated in nuclear matter. This allows 
one to examine commonly made local density approximations. Finally, we also show two 
"nonrelativistic" calculations; one with and the other without RPA correlations. Here, a 
Fermi gas is considered with an effective mass equal to the free nucleon mass. In Table ID I 
we list both sets of parameters used in the calculations. 

Ten calculations lead to somewhat complicated results. However, there are important 
motivations to examine RPA, vacuum polarization and full finite nucleus corrections to 
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simple local density approximation results. First, relativistic effects for magnetic moments, 
seen in impulse approximations, vanish in full RPA calculations Do the relativistic 

effects for A T claimed in Ref. || also vanish in RPA? Below, we see that they do not. 

Next, vacuum polarization may have an important effect on the unpolarized longitudinal 
response 
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nucleon-antinucleon pairs out of the vacuum. Therefore, the vacuum becomes better at 
screening charge. This change in vacuum polarization with M* reduces the charge visible 
to the electron probe. Relativistic RPA calculations including vacuum polarization have 
longitudinal responses reduced by about twenty five percent compared to impulse approxi- 
mations. This reduction is in good agreement with experiment for medium and heavy nuclei 
(although the experimental errors are large). In contrast, nonrelativistic calculations have 
integrated responses constrained by the Coulomb sum rule and over predict the data. Below 
we examine the effect of vacuum polarization on parity violation. 

Finally, many quasielastic calculations are done in a local density approximation (LDA). 
This is expected to be qualitatively correct for the high excitation energies of quasifree 
scattering. However, we check to see if parity violation is sensitive to any of the details of 
the response not reproduced by a LDA. In particular the high and low excitation energy 
parts of the response are poorly described in LDA. 

Figure [I] shows transverse response functions calculated in nuclear matter compared 
to data from 12 C at a momentum transfer of q = 550 MeV/c. We consider an average 
density of 12 C which corresponds to a Fermi momentum of k F = 225 MeV/c. First, free 
Fermi gas results are shown with M* = M. These are equivalent to the results reported in 
Ref. ||. Next, relativistic MFT results are shown, as in Ref. M, with an effective mass of 
M* = 0.68M. This is the self-consistent value for the MFT at a density corresponding to 
k F = 225 MeV/c. The primary effect of M* is to shift the strength to higher energies in 
agreement with data. Of course there is still substantial experimental strength at medium 
to high excitation energies. Presumably this corresponds to delta and pion production, 
multinucleon knockout and or meson exchange currents. These effects are outside the scope 
of the one-body mechanism assumed for all of our calculations. We should mention that no 
attempt was made to fit the data by selecting an optimal value for k F . As we will show 
below finite-nucleus calculations give a more realistic distribution of transverse strength. 

The contribution of other mechanisms to the quasifree response is an important source 
of uncertainty for parity violating experiments. At this time it is not well known the relative 
importance of different processes much less their isoscalar vs. isovector contribution. One 
might argue that Delta or pion exchange currents are isovector. However, there could be 
isoscalar contributions from other mesons or from two-nucleon correlations. Any effect 
from these other contributions will only strengthen our conclusion that there are important 
nuclear structure uncertainties in parity violating quasielastic electron scattering. 

We also show RPA responses in Fig. [1]. These use a simple isoscalar interaction consisting 
of sigma and omega meson exchange. The parameters are taken from the MFT fit to 
the saturation density and binding energy of nuclear matter |T2| , |TB|j . We assume that the 
isovector interaction is zero. Consequently, since the total response is dominated by the 
isovector part, RPA has almost no effect on the unpolarized transverse response. 

Next, in Fig. ^|we show the isoscalar contribution to the transverse response. In principle, 
this response contributes to the observable unpolarized cross section. However, because it 
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is so small (compare the scales of Figs. [I] and 0) it is essentially unconstrained by data. As 
M* decreases the isoscalar response is predicted to be significantly enhanced. This effect 
was discussed in Ref. |5j. 

The RPA correlations are seen to shift strength downward (due to an attractive particle- 
hole interaction) and somewhat reduce the area. This is true in both the MFT and M* = M 
calculations. Note, there is no simple relation between the MFT RPA and the free Fermi-gas 
results. This is in contrast to the situation for (isoscalar) magnetic moments where there 
was a Fermi liquid theory theorem constraining the results f20| , |2T|| . Instead, the quasifree 
response at finite momentum transfer is a dynamical quantity and depends on the details of 
the model. There is no theorem which constrains the finite-g response. 

The ratio TZ T of isoscalar to isovector transverse responses is shown in Fig. |||. This is 
the only nuclear structure information needed to determine A T (see Eq. (|37|)). In general 
there is a large enhancement in this ratio for the MFT calculations. However, there are also 
large effects from RPA correlations in both the MFT and M* = M calculations. Figure |3| 
also shows A T deduced from 1Z T by using Eq. (p7|). This quantity is directly related to 
the experimental asymmetry. The 5 to 10 percent spread in A T is much larger than the 
perhaps one to two percent accuracy needed to measure strange quark effects or radiative 
corrections to the axial contributions. 

We now consider full finite nucleus calculations for 12 C. Our finite nucleus calculations 
have been described in Ref. [|16|] and are based on calculating the nucleon Green's function 
directly in the finite system. The RPA integral equations are then solved by matrix inversion 
in momentum space. The interaction parameters are the same ones used to fit ground-state 
properties so the calculations explicitly conserve current. 

Figure ^ shows MFT results for the total transverse response. First, the finite nucleus 
calculations are qualitatively similar to the nuclear matter results. However, the shift in the 
position of the peak and the broadening of the width are not as pronounced as in the nuclear 
matter case. In addition, there is a high energy tail in the finite nucleus response which 
is not reproduced by the nuclear matter calculations. Figure |5| shows isoscalar responses. 
Again, the finite nucleus responses have a high energy tail not present in nuclear matter. 
However, the LDA results are qualitatively similar for the lower energies. 

The isospin ratio TZ T is shown in Fig. || This is qualitatively similar in finite nucleus and 
LDA calculations. However, the uncorrelated finite nucleus calculations predict an enhance- 
ment in TZ T at high excitation energies. This is due to the different energy dependences of 
n" (which is dominant for the isoscalar response) and n" (which dominates the isovector) 
in the finite nucleus results. Note, we have not yet investigated the origin of this difference. 
The A T deduced from these ratios are also shown in Fig. ^. Again there is a 5 to 10 percent 
dispersion in the results. 

Finally, we consider RHA calculations which include vacuum polarization |16|J2^] . Both 
nuclear matter and finite nucleus RHA results are shown in Figs. ([7], |8|, ^). The most 
important difference between these curves and MFT results may not be intrinsic to vacuum 
polarization. Instead, the RHA has a somewhat larger effective mass. This is because the 
RHA includes a term describing the change in the energy density of the vacuum (with M* ) . 
This term tends to resist changes in the nucleon mass. [Also, the RHA reproduces nuclear 
matter saturation with somewhat smaller meson couplings]. Therefore, RHA results tend 
to be intermediate between MFT and M* = M calculations. In particular, the isoscalar 
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response is not enhanced as much in the RHA as it is in the MFT. 

In summary we have presented results for A T for ten different calculations in Figs. (|| 
|6] and |9[). Depending on the approximation uncorrelated vs. RPA or MFT vs. RHA one 
obtains substantially different results. The five to ten percent dispersion in A T is larger than 
the accuracy needed to extract strange quark information. 

For completeness we also show axial response functions and the A T > derived from them 
in Fig. |TD|. Since we have ignored strange-quark contributions, the axial response is purely 
isovector. Thus, (isoscalar) RPA correlations have no effect on the response and we only 
report uncorrelated results. 

We see that the position and width of the quasielastic peak vary considerably reflecting 
the different effective masses of the various models. In addition, the finite nucleus calcula- 
tions show their characteristic high-energy tails. Again, we observe differences of the order 
of 5 to 10 percent in the results. 

Finally, in Figs. ( |TT| , |T2"| , and[T3|) we simply plot all ten predictions for the total asymmetry 
A = A L +A T +A T > . These show the total dispersion in the different results. The experimental 
asymmetry is simply a r times these curves. At q — 550 MeV/c and at an excitation energy 
of to = 150 MeV the overall magnitude is a r = —2.52 x 10~ 5 (see also Table |TV| ). 

IV. SUMMARY AND CONCLUSIONS 

In principle, a parity violation measurement in quasielastic electron scattering can pro- 
vide information on strange quark contributions and radiative corrections to weak axial 
currents. This could greatly enhance the value of an elastic electron-proton experiment. 
However, in order for this to be useful one must understand nuclear structure corrections to 
the quasielastic measurement. 

The quasielastic measurement is primarily sensitive to the ratio of isoscalar to isovector 
transverse responses 7Z T . We have calculated these responses in a variety of relativistic mean 
field and RPA models for 12 C and nuclear matter at a momentum transfer of q = 550 MeV/c. 
We find that 1Z T is sensitive to relativistic effects from a strong scalar field (small M*) and 
from RPA correlations. Furthermore, full finite nucleus calculations predict large changes 
in 1Z T in the high energy tail of the response. 

These nuclear structure effects could change the asymmetry A by as much as ten per- 
cent. This uncertainty is so large that a quasielastic measurement may not provide useful 
information on strange quark contributions or radiative corrections (unless these effects can 
be significantly constrained by future work). We note that a quasielastic measurement from 
deuterium rather than from a heavier nucleus should not be plagued by these uncertainties. 
We expect relativistic corrections to be much smaller for deuterium. 

Alternatively, a quasielastic parity violating measurement could provide information on 
relativistic and nuclear structure effects on the isoscalar transverse response. This small 
response is interesting and not well determined by other data. This information might be 
obtained by comparing a deuterium measurement to that from a heavier nucleus such as 
12 C. 
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NUCLEON FORM FACTOR 

We adopt the form factor parameterization used in Ref. ||. First the electromagnetic 
form factors are written in terms of a simple dipole, 

G = (1 + 4.97r)" 2 , (39) 

F™ = [1 + t(1 + \ p )]G/(1 + t), (40) 

FW = V?/(l + r), (41) 

F^= T X n (l-v)G/(l + T), (42) 

F^ = X n (l+rrj)G/(l+r). (43) 

Here the anomalous moments are, 

X p = 1.793 , A n = -1.913 , (44) 

and 

7] = (1 + 5.6T)- 1 . (45) 

This parameterization is good for the neutron form factors provided r < 1. Finally, the 
isovector axial form factor is (g a =1.26) 

G { a = —(1 + 3.53r)~ 2 . (46) 
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FIGURES 

FIG. 1. Transverse response as a function of energy loss for 12 C at q = 550 MeV. The solid 
(dot-dashed) line gives the uncorrected (RPA) Fermi gas result. The dashed (dotted) line gives the 
uncorrelated (RPA) result using an effective mass of M* = 0.68. All calculations ignore vacuum 



polarization and used a Fermi momentum of k F = 225 MeV. Experimental data is from Ref. 24]. 

FIG. 2. Same as Fig. [l] but for the isoscalar contribution to the transverse response. 

FIG. 3. Same as Fig. || but for: (a) the transverse isoscalar to isovector ratio and (b) the 
transverse contribution to the parity- violating asymmetry at 9 = 150°. 

FIG. 4. Transverse response as a function of energy loss for 12 C at q = 550 MeV. The solid 
(dot-dashed) line gives the uncorrelated (RPA) finite-nucleus result. The dashed (dotted) line 
gives the uncorrelated (RPA) result using an effective mass of M* = 0.68. All calculations ignore 
vacuum polarization and (the nuclear matter ones) used a Fermi momentum of k F = 225 MeV. 



Experimental data is from Ref. [24]. 



FIG. 5. Same as Fig. |4] but for the isoscalar contribution to the transverse response. 

FIG. 6. Same as Fig. |] but for: (a) the transverse isoscalar to isovector ratio and (b) the 
transverse contribution to the parity- violating asymmetry at 9 = 150°. 

FIG. 7. Transverse response as a function of energy loss for 12 C at q = 550 MeV. The solid 
(dot-dashed) line gives the uncorrelated (RPA) finite-nucleus result. The dashed (dotted) line 
gives the uncorrelated (RPA) result using an effective mass of M* = 0.80. All calculations include 
vacuum polarization and (the nuclear matter ones) used a Fermi momentum of k F = 225 MeV. 



Experimental data is from Ref. [24]. 



FIG. 8. Same as Fig. [7| but for the isoscalar contribution to the transverse response. 

FIG. 9. Same as Fig. ^ but for: (a) the transverse isoscalar to isovector ratio and (b) the 
transverse contribution to the parity- violating asymmetry at 9 = 150°. 

FIG. 10. Axial isovector response (a), and its contribution to the parity- violating asymmetry 
at 9 = 150° (b), as a function of energy loss for 12 C at q = 550 MeV. The solid (dot-dashed) 
line gives the MFT (RHA) finite-nucleus result. The dashed (dotted) line gives the MFT (RHA) 
nuclear matter result. Also shown are free Fermi gas results. Nuclear matter calculations used a 
Fermi momentum of k F = 225 MeV. 
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FIG. 11. Parity-violating asymmetry as a function of energy loss for C at 9 = 150° and 
q = 550 MeV. The solid (dot-dashed) line gives the uncorrelated (RPA) Fermi gas result. The 
dashed (dotted) line gives the uncorrelated (RPA) result using an effective mass of M* = 0.68. All 
calculations ignore vacuum polarization and used a Fermi momentum of k F = 225 MeV. 

FIG. 12. Parity-violating asymmetry as a function of energy loss for 12 C at = 150° and 
q = 550 MeV. The solid (dot-dashed) line gives the uncorrelated (RPA) finite-nucleus result. The 
dashed (dotted) line gives the uncorrelated (RPA) result using an effective mass of M* = 0.68. All 
calculations ignore vacuum polarization and (the nuclear matter ones) used a Fermi momentum of 
k F = 225 MeV. 

FIG. 13. Parity-violating asymmetry as a function of energy loss for 12 C at 6 = 150° and 
q = 550 MeV. The solid (dot-dashed) line gives the uncorrelated (RPA) finite-nucleus result. The 
dashed (dotted) line gives the uncorrelated (RPA) result using an effective mass of M* = 0.80. All 
calculations include vacuum polarization and (the nuclear matter ones) used a Fermi momentum 
of k F = 225 MeV. 
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TABLES 



TABLE I. Electro-weak charges for electrons and u, d, and, s quarks. Numerical values were 
evaluated using sin ~ 0.227. 



Fermion 


e(EM) g v (Weak Vector) 


<7 A (Weak Axial) 


e 


-1 -1 + 4sin 2 6\ v ~ -0.091 


+ 1 


u 


+2/3 +1 - | sin 2 6 W ~ +0.394 


-1 


d 


-1/3 -1 + | sin 2 0w ~ -0.697 


+ 1 


s 


-1/3 -1 + | sin 2 0^ ~ -0.697 


+ 1 


TABLE II. 


Isoscalar, isovector, and strange (or singlet) weak-vector and weak-axial couplings. 


Numerical values were evaluated using sin 2 8 W ~ 0.227. 




Coupling 


(Weak Vector) 


(Weak Axial) 




-4 sin 2 6 W ~ -0.909 







2-4sin 2 9 W ~ +1.091 


-2 




-1 


+1 



TABLE III. Coupling constants, meson masses, and effective nucleon mass in a mean-field 
approximation to the Walecka model. The effective nucleon mass was calculated using a Fermi 
momentum of k F = 225 MeV. 



Model g 2 s 


9l 


m s (MeV) 


m v (MeV) 


M*/M 


MFT 109.626 
RHA 54.289 


190.431 
102.770 


520 

458 


783 
783 


0.68 
0.80 


TABLE IV. Ratio of isoscalar to isovector transverse responses and parity-violating asymme- 
tries for 12 C using various nuclear-structure models at q = 550 MeV, lo = 150 MeV, and 9 = 150°. 
The nuclear matter calculations used a Fermi momentum of k F = 225 MeV. Within each model, 
the first (second) row of numbers give the uncorrelated (RPA) results. 


Model 


TZ T 


A T /a r 


A T ' 1 a a T 


A/a r 


Free Fermi Gas 


0.044 
0.030 


0.973 
1.000 


0.198 
0.201 


1.168 
1.201 


MFT[M* = 0.68] 


0.083 
0.071 


0.915 
0.937 


0.180 
0.182 


1.090 
1.117 


MFT [Fin. Nuc] 


0.072 
0.068 


0.928 
0.937 


0.184 
0.185 


1.107 
1.119 


RHA[M* = 0.80] 


0.064 
0.039 


0.943 
0.994 


0.188 
0.194 


1.126 
1.191 


RHA [Fin. Nuc] 


0.054 
0.051 


0.957 
0.969 


0.192 
0.194 


1.144 
1.164 
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